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Abstract 

We have found that the conventional generic dynamical model for gamma-ray 
bursts (GRBs) cannot reproduce the Sedov solution in the non-relativistic limit. 
Based on our refined generic dynamical model, we investigate afterglows from jetted 
GRB remnants numerically. Many new results are reached. For example, we find 
no obvious break in the optical light curve during the relativistic phase itself. But 
an obvious break does exist at the transition from the relativistic phase to the non- 
relativistic phase, which typically occurs at time 10 to 30 days. It is very interesting 
that the break is affected by many parameters, such as the electron energy fraction, 
^e, the magnetic energy fraction, the initial opening angle of the iet,6o, and the 
medium density, n. We also find that afterglows from jetted GRB remnants are 
uniformly characterized by a quick decay during the non-relativistic phase, with 
power law timing index greater than 2.1 . Afterglows from GRB 970228, 980326, 
980519, 990123, 990510 and 991208 can be satisfactorily fitted if the corresponding 
GRB ejecta are highly collimated. 



1 Introduction 

The discovery of afterglows from some gamma-ray bursts (GRBs) has opened up a new era 
in the field (Costa et al. 1997; Galama et al. 1999). The so called fireball model (Meszaros 
& Rees 1992; Sari et al. 1996), which can successfully explain the major features of GRB 
afterglows (Vietri 1997; Waxman 1997; Tavani 1997; Wijers et al. 1997), becomes the 
most popular model (see Piran 1999 and van Paradijs et al. 2000 for recent reviews). Some 
GRBs localized by BeppoSAX satellite have implied isotropic energy release of more than 
Iq5a gj^gg (^Madsen et al. 1999; Harrison et al. 1999), this has forced many theorists to 
deduce that GRB radiation must be highly collimated (Groot et al. 1998; Owens et al. 
1998; Halpern et al. 1999; Castro-Tirado et al. 1999; Sari et al. 1999; Wei & Lu 1999a, 
b). 

To tell a jet from an isotropic fireball, it seems that we must resort to the after- 
glow light curves. When the bulk Lorentz factor of a jet drops to 7 < 1/6, with 9 the 
half opening angle, the edge of the jet becomes visible, the light curve will steepen by 
^-3/4_ rpj^jg j.g^]jg^ |-]^g edge effect (Meszaros & Rees 1999; Panaitescu & Meszaros 1999; 
Kulkarni et al. 1999). Additionally, the lateral expansion of a relativistic jet will make 
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Figure 1: Velocity vs. radius for an isotropic adiabatic fireball. The 
dashed line is the familiar Sedov solution in the Newtonian phase. 
The dash-dotted line is drawn according to Eq. (1), which differs 
from the dashed hne markedly. The solid line corresponds to our 
refined model (i.e., Eq. (2)), which is consistent with the Sedov 
solution. 



the break even more precipitous (Rhoads 1999; Lamb 2000). So it is generally believed 
that afterglows from jetted GRBs are characterized by an obvious break in the light curve 
at the relativistic stage. 

In this article, I use our refined dynamical model to study the jet effect on the 
afterglow light curves. 



2 Our Refined Generic Model 

The importance of the non-relativistic phase of fireball expansion has been stressed by 
Huang et al. (1998a, b; also see Wijers et al. 1997). In the literature, it is generally 
beheved that the following equation can depict the evolution of GRB remnants (Chiang 
& Dermer 1999) 

^ = -^ (1) 

dm M ' ^ ' 

where m is the rest mass of the swept-up medium, M is the total mass in the co-moving 
frame, including internal energy. However, Huang et al. (1999) pointed out that during 
the non-relativistic phase of an adiabatic expansion, Eq. (1) cannot reproduce the familiar 
Sedov solution. This is clearly shown in Fig. 1. 



Huang et al. (1999) have proposed a refined equation 

7^ — 1 



dm + em + 2{1 — e)jm^ 



(2) 



where Mej is the initial baryon mass ejected from the GRB central engine, and e is the 
radiative efficiency. For an adiabatic fireball^ e = 0; and for a highly radiative one, e — 1. 
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Figure 2: Evolution of 7. The solid line 
corresponds to a jet with "standard" pa- 
rameters. Other lines are drawn with 
only one parameter altered: the dashed 
line corresponds to ^0 = 0.1, and the 
dash-dotted line corresponds to n — 10^ 
cm-=^ (Huang et al. 2000c). 
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Figure 3: The effect of on the R-band 
light curve. 5'r is the R-band flux density. 
The solid line corresponds to a jet with 
"standard" parameters. Other lines are 
drawn with only altered. The dashed, 
dash-dotted, and dotted lines correspond 
to — 0.03, 0.4, and 1.0 respectively 
(Huang et al. 2000c). 



Huang et al. (1999) have shown that Eq. (2) is correct for both radiative and adiabatic 
fireballs, and in both ultra-relativistic and non-relativistic phases (c.f. Fig. 1). 



3 Evolution of Jetted GRB Remnants 

Using our refined dynamical model, the evolution of the beamed ejecta can be described 
by (Huang et al. 2000a, b, c, d): 



dR 



dm 
dR 



/5c7(7 + (3) 



27ri?^(l - cos^)nmp, (4) 



dt R 

d^y 7^ — 1 



dm Mej -|- em -|- 2(1 — e)7m' 



(6) 



^^=w-l)(7-l)Y^-^c^ (7) 

where mp is the proton mass, Cg is the co-moving sound speed, 7 ~ (47 + l)/(37) is the 
adiabatic index, and c is the speed of light. Below we will consider only adiabatic jets, 
for which e = 0) . 



A strong blastwave will be generated due to the interaction of the jet and the ISM. 
Synchrotron radiation from the shock accelerated ISM electrons gives birth to afterglows. 
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Figure 4: The effect of ^| on the 
optical hght curve. The sohd hne 

corresponds to a jet with "stan- 
dard" parameters. Other hnes are 
drawn with only altered. The 
dash-dotted, dashed, and the dot- 
ted hnes are for ^| = 10-^ 10"^ 
and 10~^ respectively (Huang et 
al. 2000c). 



Figure 5: The effect of Cb on the X-ray 

light curve. Fx is 0.1 — 10 keV band flux, 
in units of ergs-cm^^-s^^. The solid line cor- 
responds to a jet with "standard" parameters. 
Other lines are drawn with only one parameter 
altered. The dash-dotted line is for — 1-3, the 
dashed line is for ^| = 0.1, and the dotted line 
is for = 1 (Huang et al. 2000d). 



As usual we assume that the magnetic energy density in the co-moving frame is a fraction 
^1 of the total thermal energy density (S'^/Stt = ^b^'), and that electrons carry a fraction 
of the proton energy. This means that the minimum Lorcntz factor of the random 
motion of electrons in the co-moving frame is 7c.mm = ^c(7 — l)"^p(p — 2)/[mc(p — 1)] + 1, 
where p is the index characterizing the power law energy distribution of electrons, and 
me is the electron mass. 



4 Numerical Results 

For convenience, let us define the following initial values or parameters as a set of "stan- 
dard" parameters: initial energy per solid angle Eq/Qq = 10^^ crgs/47r, 70 = 300, n = 1 
cm-^ ^1 = 0.01, p = 2.5, = 0.1, Oo = 0.2, ^obs = 0, L'l = lO'' kpc, where ^obs is the 
angle between the line of sight and the jet axis, and Di^ is the luminosity distance. 

We have followed the evolution of jetted GRB remnants numerically and calculated 
their afterglows (Huang et al. 2000a, c, d). Fig. 2 shows the evolution of 7 for some 
exemplary jets. We see that the ejecta will cease to be highly relativistic at time t ~ 10^ 
— 10^ s. This gives strong support to our previous argument that we should be careful 
in discussing the fireball evolution under the simple assumption of ultra-relativistic limit 
(Huang et al. 1998a, b, 1999, 2000a, b, c, d). 

Fig. 3 illustrates the effect of on the optical (R-band) light curves. In no case 
could we observe the theoretically predicted light curve steepening (with the break point 
determined by 7 ~ 1/6*) during the relativistic stage itself. The reason is: at time of 
7 ~ 1/9, the jet is already in its mildly relativistic phase and it will become non-relativistic 



soon after that, so the break due to the edge effect and the lateral expansion effect does 
not have time to emerge during the relativistic phase (Huang et al. 2000a, c). Further 
more, since 7 is no longer much larger than 1, conventional theoretical analyses (under the 
assumption of 7 ^ 1) are not proper. However, when is small, an obvious break does 
appear in the light curve, but it is clearly due to the relativistic-Newtonian transition. 
The simulation by Moderski et al. (2000) does not reveal such breaks, because their model 
is not appropriate for Newtonian expansion. When is large, the break disappears. This 
is not difficult to understand. According to the analysis in the ultra-relativistic limit, 
the time that the light curve peaks scales as oc Ce^^ld)^^^- ^ig. 3 shows this trend 
qualitatively. In the case of = 1-0, is as large as ~ 10^ s, then we can not see the 
initial power law decay (with timing index a ~ 1.1) in the relativistic phase, it is hidden 
by the peak. So the break disappears. We should also note that in all cases, light curves 
during the non-relativistic phase are characterized by quick decays, with a > 2.1. This 
is quite different from isotropic fireballs, whose light curves steepen only shghtly after 
entering the Newtonian phase (i.e., a ~ 1.3). 

Fig. 4 illustrates the effect of on the optical light curves. Interestingly but not 
surprisingly, we see that ^| has an effect similar to ^e- for small ^| values, there are 
obvious breaks at the relativistic-Newtonian transition points; but for large values, the 
break disappears, we could only observe a single steep line with a > 2.1. This is also due 
to the dependence of tm on (Huang et al. 2000c). 

We have also investigated the effects of other parameters such as Oq, n, p on the optical 
light curves (Huang et al. 2000c). When 6q > 0.4, the light curve becomes very similar 
to that of an isotropic fireball and no steep break exists. In the case of a dense medium 
(n > 10^ cm-^), the expansion becomes non-relativistic quickly, so that the break also 
disappears and we could only observe a quick decay with a > 2.1. But generally speaking, 
p does not affect the break notably. 

In Fig. 5, we show the effects of ^e^Ci 9q on the X-ray afterglow light curves. 
Again we see that a large value of any of these parameters will make the relativistic- 
Newtonian break disappear. X-ray afterglows from GRB 980519 decayed very quickly, 
with a ~ 1.83 ± 0.3. We suggest that they may be produced by a jet with a large or 
a large (Huang et al. 2000d). 

5 Comparison with Observations 

Optical afterglows from GRB 990123, 990510 are characterized by an obvious break in the 
light curve, and afterglows from GRB 970228, 980326, 980519, 991208 faded rapidly. We 
suggest that these phenomena are due to beaming effects. We have fitted these afterglows 
based on our refined jet model and find that the observations can be reproduced easily with 
a universal initial half opening angle ^0 ~ 0.1 (Fig. 6). The obvious light curve break in 
GRB 990123 is due to the relativistic-Newtonian transition of the beamed ejecta, and the 
rapidly fading afterglows come from synchrotron emissions during the mildly relativistic 
and non-relativistic phases. We thus strongly suggest that the rapid fading of afterglows 
currently observed in some GRBs is evidence for beaming in these cases. 
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Figure 6: Optical afterglows from GRB 980519, 990123, 990510, 
991208 and our best fit to them by employing the refined jet model. 
Observed data are collected from the literature. Insets show the 
evolution of 7 in our models (Huang et al. 2000b). 



6 Discussion and Conclusions 

To conclude, in this article, we try to answer the following question: does the jet effect lead 
to any breaks in the afterglow light curve? We have shown clearly that the theoretically 
predicted light curve steepen due to the edge effect and the lateral expansion effect in 
fact does not exist, consistent with recent numerical results by some other researchers 
(Moderski et al. 2000). Our new finding is that a striking break does appear in the 
light curve at the relativistic- Newtonian transition point. However, this break is affected 
by many parameters, such as ^c^b^^o and also n. Increase of any of them to a large 
enough value will make the break disappear. So, although a sharp break in the light 
curve is a good indicator for jet (note that this break is due to the relativistic-Newtonian 
transition, as pointed out above), the jet effect does not definitely lead to sharp breaks 
in the afterglow light curves. This seems to make the task of distinguishing jets from 
isotropic fireballs even more difficult. 

Fortunately we have another helpful tool. Although whether the break appears or not 
depends on parameters, afterglows from jetted GRB remnants are uniformly characterized 
by a quick decay during the non-relativistic phase, with a > 2.1. This is quite different 
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from isotropic fireballs. We strongly suggest that it can be regarded as a fundamental 
characteristic of jet and may be used to judge the degree of beaming observationally. For 
example, the rapid fading of optical afterglows from GRB 970228, 980326, 980519, 990123, 
990510 and 991208 has been argued as evidence for beaming in these GRBs (Huang et al. 
2000b). 

Showing only a single fiat line (with a ~ 1.1) in the optical light curve, radiation 
from GRB 971214 and GRB 980703 is not likely be highly coUimated. Since they have 
indicated isotropic energies of ~ 0.17 MqC^ and ~ 0.06 MqC^ respectively (Kulkarni et 
al. 1998; Bloom et al. 1998), which could hardly be met by any kind of compact stellar 
objects, we see that the energy crisis is really a problem. 
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